ffftC  Fiigcofiy 


P  *r*  ^  // 

"  '  ’  V:  ?JL' 


Cornell  University 


Magnetomechanics 

Research 


Thi,  riorum^,  har~boon~approved_ 

"  1  P  "  J I .  -  ri:loO'U'  C”H  ~n!n-  it 
distribution  id  unlimilod.  S 


Departments  of 
Theoretical  &  Applied 
Mechanics 
and 

Structural  Engineering 

8l  8  26  054 


EDDY  CURRENT  FLOWS 
AROUND  CRACKS  IN  THIN  PLATES 
FOR  NONDESTRUCTIVE  TESTING(l) 

S.  Mukherjee(2) 

M.A.  Morjaria^ 
and 

F.C.  Moon^ 

Department  of  Theoretical  and  Applied  Mechanics 
Cornell  University 
Ithaca,  NY  14850 

March  1981 


(1) 

(2) 

(3) 

(4) 


Sponsored  in  part  by  the  U.S,  Office  of  Naval  Research  under  Contract 
No.  N00014-79-C-0224. 


Associate  Professor 
Research  Associate 
Professor 


I  lus  i  in  ,  hi  fi.i.;  Ijotin  approved 
*'  1  1  :|1'  ''  '  '■  ''  asa  and  salo;  its 
I  ■  i  ,  ;  i:nlii)it{i‘d. 


SECURITY  CL  ASSIFICAT1QN  OR  THIS  RAGE  W>on  Dota  gnorodj  \  /  CT*'  '  y  I  I  | 

REPORT  DOCUMENTATION  PAGE  bmSSMSplStiWorm 

T  REPOBT NUMBEB  ]J,  GOVT  ACCF.SSION  NO,  J.  BEatPYEHT’S  CATALOG  numbeb 

AO-  flicm&n  r<?  ) 


BEPOBT  NUMBER 


i  :  EDDY  CURRENT  jXOWS  AROUND  CLACKS  IN  THIN 
I O  J*LATES  FOR  NOtTOESTFUjfcT IVE  TESTING  ,  ' 


IOO  COVERED 


Topical  report* 
Aug  |M8 0  Mara 


1 7 .  AUTHORfiJ 


y.  v  _ _ _ _ _ 

» ,4  S.  iMukherjcie^  M.A.  /Mor  jaria  land 
L  !  J|.C./Moon  M . 1 . 


It.  PERFORMING  ORGANIZATION  NAME  AND  ADDMII 


)  /'  Departments  of  Structural  Engineering  and 
J  )'  Theoretical  and  Applied  Mechanics,  Cornell 
University,  Ithaca,  NY  1A853 


L4-79-C-/224V 


PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  t  WORK  UNIT  NUMBERS 


NR  064-621 


M.  CONTROLLING  OFFICE  NAME  AND  AQORESS  y~ - '1  H  SIS 

Office  of  Naval  Research  Resident  Repr^-  /  ill  Mar«haM81  l 

sentative,  Room  323,  Federal  Bldg.,  100 - numscr  o**a<w 

_ State  Street.  Rochester.  NY  14614 _ — 8 - 

MONITORING  AGENCY  NAME  A  AOOREUCU  dltlerant  Item  Canttalltnl  Olttai)  IB.  SECURITY  CLAIt.  (ot  fhli  ttpott) 

Director,  Structural  Mechanics  Program,  TT  1  . 

Material  Science  Division,  Office  of  Naval  unciassirieu 

Research,  Arlington,  VA  22217  is*,  decl a*bii>i cation/ oownorai 


Unclassified 


A.  DECL  ASSIFICATION/  DOWN  GRADING 
SCHEDULE 


MS.  DISTRIBUTION  STATEMENT  Co  I  Ihl ,  Rtport) 


No  restriction 


I  ilx!  i1.'  I  l-'l  I  t upprovfd  j 

I  I  t,  i'll  -  ;.  I .  '  i  M  ..-I  VI  |-  ,  .1:  j 

• "  .i .; 1  ■  ii  i .  . 1 1 1 1 1 • , » i i . i i-i .  I 


MT.  DISTRIBUTION  STATEMENT  (a!  Iha  aka  treat  entered  In  Bleak  30,  II  dl  Iterant  Item  Report; 


IS.  supplementary  notes 


It.  KEY  WORDS  ('Confirm*  on  rororo*  otdo  it  nocoooory  *nd  Identity  by  block  nueifcorj 

Eddy  currents,  boundary  element  method,  nondestructive  testing, 
magnetomechanics,  numerical  methods 


20.  ABSTRACT  (Continue  on  rovouo  old*  It  nooooooty  ond  Identity  by  klooJt  numkor) 

*  The  boundary  element  method  is  used  to  calculate  the  Induced  electric 
currant  flow  around  cracks  in  thin  conducting  plates.  A  low  frequency  approx¬ 
imation  leads  to  a  Poisson  equation  for  the  current  density  potential  or 
stream  function.  A  kernel  is  used  which  produces  the  correct  singularity  at 
the  crack  tip.  The  boundary  condition  on  the  crack,  derived  from  Faraday's 
law,  requires  the  line  integral  of  the  current  density  around  the  crack  to  be 
zero.  Numerical  results  for  induced  currents  due  to  a  circular  induction  coil 

DD  ,  jan*7,  1473  edition  OB  i  Nov  «•  it  OBSOLETE  Unclassified  ' 

5/ N  0 1  02-  LF-  01 4-  6A0  I  nnniTv  CLAturiCATiau  OB  Tuii  BASE  TVhu  Dili  WntarU. 


SECURITY  CLASSIFICATION  OB  THIS  RAO!  CWNan  Gala  InlgfWj 


/o3  n  I 


apwm  Muktf cuhtne  i »  .***.’ »m#p 


SICUFITV  CLASSIFICATION  OF  THIS  FAO*  fWN««  Dm*  tnlmn*) 


SUCUNITV  CLASSIFICATION  OF  THIS  FAC«f*IH«l  B«* 


INTRODUCTION 


Induced  electric  currents  are  generated  in  conductors  by  time  vary¬ 
ing  magnetic  fields.  When  the  source  of  the  field  is  outside  the  body, 
the  induced  currents  must  flow  in  closed  paths,  hence  the  designation 
"eddy  currents".  In  most  problems  eddy  currents  are  unwanted  since  they 
are  a  source  of  heat  and  energy  loss  and  in  some  applications  can  cre¬ 
ate  dynamic  forces  and  magnetic  pressures.  A  few  applications,  such  as 
magnetic  forming  and  levitation,  have  exploited  the  dynamic  force  pro¬ 
ducing  capability  of  eddy  currents.  A  third  interest  in  eddy  currents  is 
their  potential  for  nondestructive  testing.  The  presence  of  flaws  or 
cracks  interrupt  the  natural  flow  of  electric  current,  and  the  detection 
of  the  change  of  electron  flow  can  give  a  clue  to  the  presence  of  flaws 
in  solid  conductors. 

The  calculation  of  eddy  currents  in  conductors  is  generally  carried 
out  by  either  using  a  magnetic  or  an  electric  potential.  A  comparison 
of  the  two  methods  is  given  by  Carpenter  [l].  The  electric  field  or 
current  density  potential  has  the  advantage  that  it  need  be  calculated 
only  in  the  conductor,  whereas  the  magnetic  potential  must  be  solved  both 
inside  and  outside  the  conductor.  The  latter  method,  therefore,  poses 
potential  problems  for  numerical  methods. 

Numerical  methods  must  generally  be  used  for  the  solution  of  eddy 
current  problems  in  conductors  of  complex  shape.  The  finite  element  meth¬ 
od  (.TEH)  and  the  discrete  circuit  element  method  have  been  used  for  some 
years  for  the  solution  of  these  problems.  Recently,  the  boundary  element 
method  (BEM)  (also  called  the  boundary  integral  equation  method)  lias 
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been  applied  to  problems  in  electromagnetics.  Wu  et  al  [2]  and  Ancelle 
et  al  [3]  have  addressed  magnetostatic  problems  by  the  BEM  while  Trowbridge 
[It]  has  considered  magnetostatic  problems  and  eddy  current  problems  by 
the  magnetic  potential  approach.  Very  recently,  Salon  and  Schneider  [5] 
have  solved  problems  of  eddy  current  flow  in  long  prismatic  conductors 
by  the  BEM  based  on  an  electric  potential  approach.  The  boundary  ele¬ 
ment  method  has  the  important  advantage  that  only  the  boundary  of  a  body 
(rather  than  the  entire  domain)  needs  to  be  discretized  in  a  numerical 
solution  procedure  -  thus  effectively  reducing  the  dimension  of  a  problem 
by  one.  However,  a  full  matrix  must  be  treated  in  the  BEM  while  the 
FEM  requires  operations  on  sparse  matrices. 

The  direct  boundary  element  approach  [2-5]  uses  a  singular  solution 
of  a  differential  equation  in  an  infinite  domain  as  a  kernel  in  the  cor¬ 
responding  integral  equation.  This  direct  approach  can  be  used  in  simply 
connected  as  well  as  multiply-connected  domains.  However,  if  a  cutout  in 
a  conducting  plate  is  a  crack,  numerical  difficulties  might  arise  from 
discrete  modelling  of  the  crack  boundary,  This  difficulty  can  be  over¬ 
come  if  modified  kernels  are  used  so  that  the  new  kernels  are  the  singu¬ 
lar  solutions  of  the  governing  differential  equations  for  an  infinite  re¬ 
gion  with  a  crack  already  present  in  it.  This  technique  has  been  recent¬ 
ly  developed  for  two-dimensional  harmonic  and  biharmonic  operators  in 
connection  with  study  of  stresses  near  crack  tips  in  bodies  undergoing 
inelastic  deformation  [6-8].  Use  of  these  modified  kernels  allows  the 
proper  boundary  conditions  to  be  satisfied  exactly  over  the  entire  crack 
surface  and  discretization  of  the  crack  surface  is  no  longer  necessary 
in  a  numerical  solution  procedure.  The  method  is  thus  perfectly  suited 
to  the  study  of  two-dimensional  problems  of  eddy  current  flow  in  cracked 


bodies. 
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The  purpose  of  this  paper  is  the  study  of  eddy  currents  in  thin 
cracked  plates  with  a  view  towards  detection  of  cracks  or  flaws  ty  non¬ 
destructive  testing.  An  analytical  formulation  using  an  electric  poten¬ 
tial  is  first  presented  for  the  determination  of  eddy  ourrent3  in  a  thin 
flat  conducting  plate  with  a  line  crack  present  in  it.  The  applied  mag¬ 
netic  field  is  assumed  to  be  harmonic  in  time  but  can  have  an  arbitrary 
spatial  distribution  inside  the  plate.  This  results  in  Poisson’s  equa¬ 
tion  for  the  electric  potential.  A  boundary  element  formulation  is  next 
presented  using  modified  kernels  for  the  Laplaeian  operator  in  two  dimensions. 

Numerical  results  are  given  for  eddy  currents  in  a  center-cracked 
square  plate  with  the  applied  field  being  that  due  to  a  circular  coil. 

Stream  lines  are  given  for  various  positions  of  the  coil  relative  to  the 
crack.  The  induced  temperature  at  any  point  in  the  plate  is  proportional 
to  the  square  of  the  density  of  the  induced  current  at  that  point.  Calcu¬ 
lated  induced  temperature  profiles  are  presented  for  various  coil  posi¬ 
tions.  An  eddy  current  Intensity  factor,  analogous  ^o  a  stress  intensity 
factor,  1b  defined  at  a  crack  tip.  Finally,  experimental  results  are 
presented  for  an  infrared  isotherm  of  induced  eddy  currents  in  a  cracked 
aluminum  plate. 


GOVERNING  DIFFERENTIAL  EQUATIONS 

A  thin,  flat  uniform  plate  made  of  a  conducting  material  is  shown 
in  Fig.  1.  The  plate  boundary  can  be  arbitrary,  its  thickness  is  h 
and  the  conductivity  of  the  plate  material  is  c.  The  plate  has  a  line 
crack  of  length  c  ■  2a  present  in  it.  The  crack  can  have  arbitrary 
orientation  relative  to  the  outside  boundary  of  the  plate. 


The  coordinate  system  chosen  is  shown  in  Fig,  1.  The  x^  and  x2 
axes  lie  on  the  mid-surface  of  the  plate,  with  the  x^  axis  along  the 
crack  and  x2  normal  to  it.  The  origin  of  coordinates  lies  at  the  cen¬ 
ter  of  the  crack  in  the  widsurface  of  the  plate.  The  x^  axis  is  normal 
to  the  plate.  Also,  the  unit  vectors  i,  ^  and  k  are  oriented  along 
the  coordinate  axes. 

Consider  a  current  density  J,  which  is  induced  in  the  plate  by  an 
oscillatory  magnetic  field  B°  outRide  the  plate.  The  current  distri¬ 
bution  is  assumed  to  he  uniform  across  the  plate  thickness  and  oscillatory 
in  nature.  The  skin  depth,  which  is  inversely  proportional  to  the  square 
root  of  the  frequency,  is  assumed  to  he  large  compared  to  the  plate  thick¬ 
ness.  Under  these  assumptions  no  bending  occurs  in  the  plate. 

According  to  Ohm's  law 

J  -  aE  (l) 

aw  At 

where  E  is  the  electric  field  (the  Hall  effect  or  magnetoresistive  terms 
are  neglected  in  Ohm's  law). 

For  low  frequency  currents,  the  continuity  condition  is 

7-J  -  0  (2) 

where  7  is  the  gradient  operator  in  two  dimensions.  Thus,  a  stream 
function  (or  electric  potential)  ^(x^.Xg)  can  be  defined  such  that 

J  ■  7x(^k)  ■  -kx7ti»  (3) 

so  that  J2  -  - 
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Using  Faraday's  law  of  induction 


7  x  E 


3B 

at 


(M 


with  B  the  total  magnetic  field  inside  the  plate  and  t  time,  the 
governing  differential  equation  for  the  stream  function  is 


3t 


»  £  (B!  *  B3> 


In  the  above,  is  the  self  magnetic  field  inside  the  plato  due 

to  the  current  J.  In  general,  this  field  can  be  obtained  from  the  Biot- 
Savart  law  as  an  integral,  over  the  plate,  of  a  kernel  times  the  stream 
function  ijj.  Equation  (5)  would  then  become  an  integro-differential 
equation  for  the  stream  function  [9].^/lf,  however,  the  applied  field  is 
sinusoidal  and  the  resulting  skin  depth  is  greater  than  ten  times  the  plate 
thickness,  the  self  field  term  can  be  neglected  relative  to  the  applied 
field  B°  [10],  Under  this  assumption,  and  with  B^  ■  B°  e^ut  (where 
1  ■  t£>l  and  ui  is  the  oscillation  frequency) ,  the  spatial  part  of  the 
stream  function  ip  satisfies  the  equation 


iuoB^ 


f(x1,x2) 


(6) 


which  is  a  two-dimensional  Poisson's  equation  with  a  prescribed  nonhomo- 
geneous  term.  For  simplicity,  the  same  notation  is  used  in  the  follow¬ 
ing  for  the  amplitudes  of  the  various  oscillatory  functions,  as  has  been 
used  so  far  for  the  functions  themselves. 

Eddy  current  distributions  using  a  similar  electric  potential  approach 
have  been  obtained  by  Salon  and  Schneider  [5].  Their  formulation  is 
valid  for  the  determination  of  currents  in  long  prismatic  conductors 
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and  results  in  a  Helmholtz  equation  for  the  stream  function  0.  Thus, 
the  work  of  Schneider  and  Salon  is  analogous  to  plane  strain  problems 
in  mechanics,  while  the  present  paper  addresses  conductors  in  the  shape 
of  thin  plates.  This  is  analogous  to  plane  stress. 

The  current  must  be  tangential  to  the  boundary  of  the  plate  at  a  point 
on  it.  Thus,  for  a  point  on  either  the  crack  boundary  30^  or  on  the 
outside  boundary  3Cg  (Fig.  l) 

J.n.fJ--°  {(  (T) 

where  n  is  an  unit  normal  to  the  boundary  at  a  point  on  it  and  s  is 
the  distance  measured  along  a  boundary  in  the  anticlockwise  sense.  ThuB, 
if  \|<  is  a  constant  a.^  on  3C^  and  another  constant  a2  on  3C2, 
equation  (T)  is  satisfied.  One  of  these  constants  can  be  set  to  zero 
without  loss  of  generality  and  the  other  one  is  determined  by  the  auxil¬ 
iary  condition 

$  J«tds  *0  (8) 

3C1  '  - 

where  t  ia  an  unit  tangent  to  3C^  at  a  point  on  it.  Physically, 
this  equation  implies  that  the  net  flux  following  through  the  crack  is 
zero. 

The  boundary  conditions  on  iji,  in  this  formulation,  are  therefore 
t|>  *  0  on  the  crack  boundary  3C^  (9) 

■jjj*  *  0  on  the  outside  boundary  3C2  (10) 


(11) 
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Equations  (9-ll)  together  with  the  field  equation  (6)  constitute 
a  well  posed  problem. 

It  should  be  noted  that  this  formulation  assumes  that  no  current 
flows  across  the  crack  or  crack  tip.  This  formulation  leads  to  a  cur¬ 
rent  density  singularity  at  the  crack  tip,  as  does  analogous  formulations 
for  the  stress  at  a  crack  tip.  Physically  we  suspect  that  there  is  a 
finite  resistance  or  current  leakage  across  the  crack  tip  which  would 
relieve  the  singularity  in  actual  conductors.  However  this  is  not  con¬ 
sidered  in  this  paper.  Instead  we  will  characterize  the  current  at  the 
crack  tip  by  a  current  density  intensity  factor  analogous  to  that  in 
fracture  mechanics. 


BOUNDARY  ELEMENT  FORMULATION 


Integral  equations 

An  integral  equation  formulation  for  Poisson's  equation  (6)  can  be 
written  as  (Pic.  l) 

2ir$(p)  -  j  K(p,Q)0(Q)ds 

3C2  y 

+  JAK(p,q)f(q)dA4  (12) 

This  is  a  single  layer  potential  formulation  where  0,  a  source 
strength  function  on  the  outside  boundary,  must  be  determined  from  the 
boundary  condition  on  it  (equation  10).  The  points  p  (or  P)  and 
q  (or  Q)  are  source  and  field  points,  respectively,  with  capital  let¬ 
ters  denoting  points  on  the  boundary  of  the  body  and  lower  case  letters 
denoting  points  inside  the  body.  The  area  of  the  body  B  is  denoted 
by  A. 
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The  kernel  K(p,q),  for  a  simply  connected  region,  is  normally  chosen 
to  he  a  singular  solution  of  Laplace's  equation  in  appropriate  dimensions, 

A 

eg.  K  ■  tor  ■  Het<frU,z  )1 

Jr  w 

A 

with  $(z,z  )  ■  to(z-z  )  . 

Here  r  _  is  the  distance  between  a  source  point  p  and  a  field  point 

pq 

q,  Be  denotes  the  real  part  of  the  complex  argument  and  z  and  zq 
are  the  source  and  field  points,  respectively,  in  complex  notation. 

In  this  problem,  however,  the  kernel  must  be  chosen  such  that  it 
vanishes  on  the  crack  boundary  3C^,  Thi3  is  achieved  by  augmenting 

A  » 

4>  with  a  second  piece  4  which  equals  the  negative  of  $  when  the  source 

/  ,  * 

point  z  lies  on  3C1  (Fig.  1).  Furthermore,  $  must  satisfy  Laplace's 

equation  and  be  regular  inside  the  body  B.  For  an  elliptical  cutout 

« 

30^,  b  is  derived  by  making  use  of  the  mapping  function 

z  -  n»U)  -  |  +  m;  (13) 

which  transforms  the  region  on  and  outside  an  ellipse  in  the  z  plane 
to  a  region  on  and  inside  an  unit  circle  in  the  ?  plane.  The  parameter 
m  equals  (a-b)/(a+b)  (with  (a+b)  *  2)  in  terms  of  the  semi-major 
and  minor  axles,  a  and  b  respectively,  of  the  ellipse.  For  the  line 
crack  in  this  problem,  a  is  taken  to  be  equal  to  2  and  b  is  zero. 

Thus,  m  equals  1.  Using  this  value  of  m,  the  augmented  function 
4  is  determined  as  [6] 

b(z,z,zQ)  -  ind-rj/O  -  tnd-rjT)  (lM 
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z  ±  Jz'-\ 

where  ^  *  — — £ — 2 —  ,  ]  ri  |  <_  1 


and  the  kernel  K  in  equation  (12)  is 


K(p,q)  ■  ReU(z,z,zQ]  (15) 

A  superposed  bar  denotes,  as  usual,  the  complex  conjugate  of  a  complex 
quantity. 

Use  of  K  from  equation  (15)  in  equation  (12)  satisfies  equation 
(9).  It  has  been  proved  in  ref.  [6]  that  this  formulation  also  satis¬ 
fies  the  integral  condition  (11)  on  the  crack  surface.  Thus,  the  proper 
boundary  conditions  on  the  crack  surface  are  satisfied  in  an  implicit 
manner  and  discretization  of  the  crank  boundary  is  not  necessary  in  a 
numerical  solution  procedure. 

The  remaining  boundary  condition  (3,0)  on  the  outside  surface  is  satis¬ 
fied  by  using  a  differentiated  version  of  f 12 )  and  taking  the  limit  as  p 
inside  B  approaches  a  point  P  on  3Cg.  Defining 

“i  ■  “*<i £>  ■  -  |f) 

2  dZ 

H 2  --»«<#■)  '  -*•(&♦**> 

1  3z 

and 

*  Re [4^]  *  H.n.  (i  summed  over  1,  2) 

as  ds  l  l 

where  n^  are  the  components  of  the  unit  outward  normal  to  5C^  at 
some  point  on  it,  the  boundary  conditon  (10)  becomes 
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0-0  H. (P,Q)n, (P)G(Q)dsQ 
3C2  h 

+  /AHi(P,q)n1(P)f(a)dAq  (17) 

Equation  (17)  is  valid  for  a  point  P  on  3Cg  where  it  ia  local¬ 
ly  smooth.  It  can  he  shown  that  in  this  case  there  is  no  extra  term  due 
to  a  residue  from  the  limiting  process.  The  boundary  integral  must  be 
interpreted  in  the  sense  of  a  Cauchy  principal  value. 

The  current  componenets  and  Jg  at  a  point  p  inside  the  body 

are  obtained  from  equations  like  (12)  with  the  kernel  K  replaced  by 
H^(p,Q)  and  Hg(p,Q)  respectively.  As  p  approaches  a  point  on  3Cg 
it  can  be  shown  that 

if  2irJj  (p#)  -  hj(p*) 

then  2irJj(P*)  -  hj(P*)+  »t ^  (P*)G(P*) 

ft  # 

where  p  is  infinitesimally  close  to  P  ,  3Cg  is  locally  smooth  at 

ft  # 

P  and  tj(P  )  are  the  components  of  the  unit  anticlockwise  tangent 

« 

vector  to  3C2  at  P  .  Thus,  the  residue  is  zero  for  the  normal  component 
and  equala  G(P)/2  for  the  tangential  component  of  the  current,  as  a  point 
approaches  the  outer  boundary.  The  actual  current  components,  of  course, 
remain  continuous  as  a  point  approaches  the  boundary. 

Discretization  of  equations  and  solution  strategy 

The  outer  boundary  of  the  body,  3Cg,  is  divided  into  Mg  straight 
boundary  elements  using  (N^  ■  Ng)  boundary  nodes  and  the  interior 
of  the  body,  A,  is  divided  into  n^  triangular  internal  elements. 

A  discretized  version  of  equation  (17)  is 
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where  la  the  point  P  where  it  coincides  with  a  node  M  at  a  cen¬ 
ter  of  a  boundary  segment  on  9Cg  and  As^  and  AA^  are  boundary  and 
Internal  elements  respectively. 

A  simple  numerical  scheme  is  used  in  which  the  source  strengths 
G  are  assumed  to  be  piecewise  uniform  on  each  boundary  segment  with 
their  values  to  be  determined  at  the  nodes  which  lie  at  the  centers  of 
each  segment.  The  integrals  of  on  boundary  elements  are  evaluated 
analytically  for  the  singular  and  by  Gaussian  quadrature  for  the  regu¬ 
lar  portions.  Nonsingular  area  integrals  of  known  integrands  over  tri¬ 
angular  internal  cells  are  evaluated  by  Gaussian  quadrature.  This  in¬ 
tegrand  becomes  singular  when  the  source  point  lies  on  the  side  of  or 
inside  of  a  triangle  over  which  the  integral  is  being  evaluated.  In  fact, 
the  location  of  an  internal  source  point,  in  equation  (12),  can  be  ar¬ 
bitrary.  In  any  case,  for  singular  integrands,  the  order  of  the  singularity 
is  1/r,  r  being  the  distance  between  the  source  and  field  points. 

The  area  element,  however,  is  rd9dr  (in  polar  coordinates).  This  form 
of  dA^  is  used  to  change  the  integrand  into  a  regular  function,  and 
then  the  integral  is  evaluated  appropriately  by  Gaussian  quadrature. 

Substitution  of  the  piecewise  uniform  source  strengths  into  equa¬ 
tion  (l8)  and  carrying  out  of  the  necessary  integrations  leads  to  an 
algebraic  system  of  the  type 


(0}  -  [A]{G>  +  {d} 


(19) 
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The  coefficients  of  the  matrix  [A]  contain  boundary  integrals  of 
the  kernel.  The  vector  {d}  contains  contributions  from  the  area  Integrals 
and  the  vector  (0}  the  unknown  source  strengths  at  the  boundary  nodes. 

The  dimension  of  {&}  depends  only  on  the  number  of  boundary  elements 
on  3Cg  and  the  internal  discretization  is  necessary  only  for  the  evalu¬ 
ation  of  Integrals  with  known  integrands. 

Equation  (12)  for  the  stream  function  \p  and  analogous  equations 
for  the  current  components  and  Jg  are  discretized  in  Bimilar  fashion. 

The  solution  strategy  is  as  follows.  The  matrix  [A]  and  vector 
(d)  in  equation  (19)  are  first  evaluated  by  using  the  appropriate  ex¬ 
pressions  for  the  kernels  and  the  prescribed  function  f  in  equation 
(6).  Equation  (9)  is  solved  for  the  vector  (0),  This  value  of  {0} 
is  now  used  in  a  discretized  version  of  equation  (12)  to  obtain  the  values 
of  the  'stream  function  t|i  at  any  point  p.  Finally,  the  current  ■vec¬ 
tor  at  any  point  is  obtained  from  equations  analogous  to  (12). 

NUMERICAL  RESULTS 

Field  Due  to  a  Circular  Induction  Coll.  A  center-cracked  square  plate 
with  a  circular  induction  coil  placed  above  it  is  shown  in  Fig.  2.  The 
square  plate  is  of  side  L  with  a  center  crack  of  length  c  «  2a.  The 
coordinate  system  is  shown  in  Figs.  1,  2.  The  coil  is  of  radius  aQ 
with  its  center  at  the  point  (x^,Xg,h0).  The  induced  field  B°  at 
a  point  q  (x^.Xg)  in  the  plate,  from  the  Biot-Savart  law,  is 

,o  ^o1  j  d!x? 

?  "  TtT  $coil  ~ 


(20) 
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where  u0  is  the  permeability  of  vacuum,  I  is  the  current  in  the  coil, 

ds  a  length  element  along  the  coil  (ds  «  e,a  d<f>)  and 
«*  „  ■*  o 

-  R  ■  (x^-x1+aocoa$)  i  +  (Xg-Xg+agSimjs);)  +  bQk 
with  R  ■  |R|  . 

Thus,  the  function  fCx^iXg)  in  equation  (6)  is 
ituffu  I  ck*(ds*R) 

f(xltx2)  -  c4ff  $ooil  ^3  ^ 

The  dimensionless  line  integral  1  can  be  written  as 

k* (dsxR)c  2ir  _ °^ao^xl~xI^cos^xg~x2^ain<^'l’ao^(1<fr 

$coil  g3  o  [(x^_xi)2+(x®_Xg)2+2ao[(x°-x1)coBi(i+(x°'*x2)8in()>]+a^+h^]': 

This  integral  is  evaluated  by  Gaussian  integration  in  the  numerical 
calculations  (6  Gauss  points  between  0  and  ir/2). 

Nondlmensional  Equations.  Equation  (6),  with  f(x^,x2)  defined  by  equa¬ 
tion  (2l),  can  be  nondimensionalized  to  the  form 

,*  Ap  A  A  A 

V>(x1,x2)  -  I  (22) 

where  J  -  ,  x ^  •  kx^/c  ,  x°  ■  kx°/c 

(i«l,2)  ,  L  ■  l*L/c  ,  aQ  ■  ^aQ/o  ,  hQ  ■  ltho/c 

Afl  ■!  £  2 

v  ■  +  ^  "  - S’  *nd  the  skin  depth 

3x^  3Xg  2iri“ 


Further,  the  dimensionless  current  density  is 


Geometrical  Parameters  and  Boundary  Mesh.  The  values  of  the  geometrical 


parameters,  used  in  the  numerical  calculations,  are 

A  A  A 

L  ■  20  ,  aQ  ■  4  and  1  ,  ho  ■  1  , 
x®  ■  0,  1,  2,  3,  4,  5,  and  6  ,  xg  ■  0 

The  dimensionless  crack  length  here  is  4.  The  results  for  any  crack 
length  c  can  he  determined  from  the  above  equations. 

A  typical  boundary  element  mesh  for  the  problem  is  shown  in  Pig.  3. 

Only  the  upper  half  of  the  plate  is  modelled  due  to  symmetry.  This  mesh 
has  42  boundary  segments  and  128  internal  cells.  The  boundary  ele¬ 
ments  are  uniformly  distributed  along  the  boundary  of  the  plate.  The 
symmetry  line  is  not  discretized.  The  unknown  values  of  the  Bource  strengths 
G  (equation  18)  lie  only  on  the  boundary  nodes.  The  internal  discretiza¬ 
tion  is  necessary  only  for  the  evaluation  of  the  area  integral  in  equation 
(18).  The  integrand  in  this  integral  is  completely  known.  Nonsingular 
integrals  are  evaluated  by  Gaussian  quadrature  with  6  Gauss  points 
on  a  boundary  segment  and  T  Gauss  points  in  an  internal  cell.  A  4x4 
grid  is  used  for  the  evaluation  of  singular  area  integrals  (see  below 
equation  (3,8)). 

The  uniform  field  problem  is  analogous  to  the  Saint  Venant  prob¬ 
lem  of  torsion  of  long  prismatic  bars  with  end  couples.  The  computer 
program  used  here  has  been  verified  by  solving  torsion  problems  in  the 
absence  of  cracks  [ 11 ]  and  the  logic  for  inclusion  of  the  crack  has  been 
verified  by  solving  Mode  III  crack  problems  [6]. 

The  boundary  Integral  algorithm  was  also  applied  to  a  plate  with  a 
notch  cut  instead  of  a  crack.  The  results  were  compared  with  an  analysis 
using  a  finite  element  method,  [10],  and  the  agreement  was  very  good. 

■  ..  . . . . .  .  .1.-^ -  '  u . T"  ~  ' 
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Eddy  Currant  and  Temperature  Linas.  Eddy  current  stream  lines  (constant 
ip  lines)  are  shown  in  Figs.  4-5  for  a  coil  of  radius  4.  Fig.  4  shows 
the  lines  in  a  plate  without  a  crack  in  it,  and  Fig.  5  shows  how  the  stream 
lines  are  affected  by  the  crack  for  different  coil  positions.  A  close-up 

a  a 

of  the  stream  lines  near  a  craok  tip  for  ■  2  is  shown  in  Fig.  6.  The 
crowding  of  stream  lines  near  the  crack  tip  leads  to  large  gradients  of 
(and  therefore  large  induced  currents)  in  this  region.  The  local  tem- 

A  A 

perature  is  proportional  to  the  square  of  the  current  density  (J*J) 
and  this  leads  to  a  hot  spot  at  the  crack  tip.  This  is  shown  in  Fig.  7 
which  shows  lines  of  constant  induced  temperature.  The  contour  lines 
go  off  scale  as  one  approaches  the  crack  tip.  The  behavior  of  the  singu¬ 
larity  at  the  crack  tip  is  dismissed  later  in  the  subsection  entitled 
'Eddy  Current  Intensity  Factor'. 

Finally,  the  stream  lines  for  a  smaller  coil  (of  diameter  equal 
to  half  the  crack  length)  centered  at  (2,0,1)  are  Bhovrn  in  Fig.  0, 

Temperature  Scans.  A  matter  of  considerable  interest  in  this  approach  to 
nondestructive  testing  is  the  existence  of  hot  spots  due  to  the  presence 
of  a  crack.  Figure  9  shows  calculated  temperature  profiles  along  a  line 
slightly  above  the  craok  (xg«0.05)  /or  different  coil  positions  (x^  ■ 
0,1,2  and  6).  The  coil  radius  hare  is  four  (equal  to  the  crack  length). 

Hot  spots  are  seen  near  the  crack  tips.  The  strongest  hot  spots  arise 
when  an  edge  of  the  coil  is  near  a  crack  tip.  These  temperatures  are  much 
higher  than  other  moderate  hot  spots  elsewhere  in  the  plate.  A  discus¬ 
sion  of  experimental  results  is  given  in  the  next  section. 
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Self  Induced  Field  at  Center  of  Coll.  One  method  of  electromagnetic  non¬ 
destructive  testing  uses  one  or  more  passive  sensing  coils  together  with 
an  active  induction  coil  [12].  The  purpose  of  the  sensing  coils  is  to  meas¬ 
ure  the  self  induced  field  (hack  e.m.f.)  and  to  try  to  observe  changes 
in  hack  e.m.f.  due  to  the  presence  of  cracks.  With  this  in  view,  the 
self  induced  field  was  calculated  at  the  center  of  the  induction  coil 
for  various  coil  positions.  The  method  used  is  the  Biot-3avart  Law  which 
gives  the  induced  field  as 

t  U  h  ,  (JxR) 

B3  “  1ST  plate  *  '  —  <*3) 

where  now  the  induced  current  density  in  the  plate,  must  be  used.  The 
integral  must  be  evaluated  over  the  plate  with  dA  an  area  element  in 
the  plate.  The  current  density  is  assumed  to  be  piecewise  uniform  over 
each  internal  cell  in  this  approximate  calculation,  with  the  value  de- 
termined  at  the  centroid  (x^,Xg)  of  the  cell.  Thus,  in  this  case 

R  "  +  (x°-x°)J  +  hQk 


„I  _  uoh  r1  <iAi^Jl^x2~x2^"J2^xl“xl^ 
B3  “  T7  A  3 


where  dA^  is  the  area  of  the  ith  triangular  element. 

A  plot  of  normalised  with  respect  to  coil  position  is  shown  in 
Fig.  10.  The  values  are  normalized  with  respect  to  when  the  coil 
center  is  directly  above  the  crack  center  (i.e.  x£  ■  x°  *  0).  It  is 
seen  that  in  this  example  the  position  of  the  coil  relative  to  the  crack 


causes  little  variation  in  the  induced  field  at  the  coil  center.  From 
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these  calculations,  it  appears  that  for  low  frequencies  the  "back  e.m.f. 
method  is  not  useful  for  detection  of  cracks  and  that  the  temperature 
scan  approach  appears  to  be  much  more  promising  for  nondestructive  test¬ 
ing  In  this  case. 

Eddy  Current  Intensity  Factor.  It  is  well  known  from  elastic  fracture 

mechanics  that  stress  components  exhibit  a  square  root  singularity  near 

a  crack  tip.  It  is  therefore  expected  that  the  components  of  the  current 

vector,  in  this  problem,  should  display  similar  behavior  near  a  crack  tip. 

*2 

This,  in  fact,  is  the  case,  as  shown  in  Fig.  11  where  J  is  plotted 

A 

as  a  function  of  r,  the  distance  from  the  right  crack  tip.  The  plot 
is  for  the  coil  center  at  (1,0,1).  The  eddy  current  density  squared 
is  seen  to  be  inversely  proportional  to  the  distance  r. 

An  eddy  ourrent  intensity  factor  M^,  analogous  to  the  stress 
intensity  factor  for  Mode  III,  is  defined  here  as 

?  ■  «,,!  *  (25) 

A  plot  of  Mp.j  at  the  two  craok  tips,  aB  functions  of  coil  posi¬ 
tion,  is  shown  in  Fig.  12.  The  eddy  current  intensity  factor  is  seen 
to  peak  when  an  edge  of  the  coil  is  near  a  oraak  tip. 

Computing  Times.  All  the  computing  reported  in  this  paper  was  carried 
out  on  an  IBM  370/168  computer  at  Cornell  University.  A  typical  compu¬ 
ting  time  for  stream  lines  in  a  cracked  plate  for  a  fixed  coil  position 
(eg.  Fig.  3  )  was  ioo  e.p.u.  secondB. 
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EXPERIMENTAL  RESULTS 


Infrared.  Experiments 

Conventional  eddy  current  nondestructive  techniques  use  a  small  in¬ 
duction  coil  and  search  coils  to  induce  eddy  currents  near  tae  surface  of 
solids  and  to  measure  to  hack  emf  generated  by  these  currents  [12].  To  detct 
a  flov  or  crack,  the  coils  must  be  moved  over  the  surface  near  the  flow 
in  order  to  measure  a  change  in  voltage  in  the  search  coil.  Recently  a  new 
method  has  been  proposed  using  infrared  scanning  technology  [10].  This  method 
is  based  on  the  fact  that  eddy  currents  create  heat  and  that  this  small 
temperature  change  can  be  detected  using  an  infrared  sensitive  device. 

In  the  present  experiments  an  aluminum  plate  15  cm  x  30  cm,  0.51  mm 
thick,  had  a  crack  placed  in  its  center,  parallel  to  the  15  cm  width.  The 
crack  was  created  by  scoring  a  line  in  the  aluminum  with  a  sharp  edge  and 
flexing  the  plate  until  fatigue  produced  a  thorough  crack  in  the  plate. 

The  length  of  the  crack  was  6  cm.  One  tip  of  the  crack  was  3  cm  from  one 

edge  of  the  plate  and  the  other  was  6  cm  from  the  other  edge. 

The  induction  coil  was  wound  from  10  turns  of  copper  wire  on  a  5.1 
cm  coil  form  so  that  the  mean  diameter  of  the  coil  was  about  5.7  cm,  slight- 
lv  smaller  than  the  crack  length.  The  width  of  the  coil  was  1.3  cm  and  the 
coil  face  was  placed  6.4  mm  from  the  plane  of  the  cracked  plate. 

Pulsed  electric  currents  of  the  order  of  9.3  KA  peak  aurrent  and  3.2 

msec,  duration  were  used  with  a  rise  time  to  peak  of  about  0.7  msec.  The 
sensitivity  of  the  infrared  system  was  0.1  -  0.2°C.  If  an  infrared  scan 
is  made  of  the  plate  immediately  after  the  firing  of  the  current  pulse, 
heat  conduction  may  be  neglected,  and  the  measured  rise  :Ln  temperature 
is  proportional  to  the  integral  of  J  over  time.  (See  e.g.  [l0].) 

The  infrared  system  used  for  these  experiments  is  a  UTI  Corp.  Spec- 
trotherm  infrared  scanning  system.  Radiation  from  different  points  in 
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the  plane  of  foous  is  detected  by  a  photoconductive  crystal.  A  two- 
dimensional  scan  is  obtained  by  two  sets  of  rotating  mirrors.  The  output 
can  be  displayed  either  on  a  grey  scale  tube  or  the  voltage  can  be  color 
quantized  into  ten  colorB.  Figure  13  shows  a  black  and  white  photograph 
of  two  color  quantized  infrared  scans.  Each  of  the  scans  represent  J2  "iso 
therms"  for  different  induction  coil  positions.  The  top  photo  corresponds 
to  a  coil  position  centered  at  the  middle  of  the  crack.  The  bottom  photo 
corresponds  to  an  induction  coil  position  centered  to  the  left  of  the  crack 
with  the  coil  center  U  cm  from  the  left  edge  of  the  plate,  These  photos 
show  that  a  "hot"  spot  forms  at  the  tip  of  the  oraok  due  to  the  current 
flowing  around  the  crack.  (The  white  line  shows  the  vertical  position 
of  the  crack.)  In  Figure  13  a  right  crack  tip  shows  up  an  circular  iso¬ 
therms.  When  the  aoil  position  is  moved  to  the  left  the  left  crack  tip 
shows  up  as  a  hot  spot,  whioh  appear  as  circular  isotherms.  The  hot  spot 
on  the  left  edge  is  due  to  increase  of  current  density  near  the  edge. 

Needless  to  say,  the  color  photographs  are  more  dramatic,  But  the  exper¬ 
iments  show  the  same  qualitative  behavior  as  the  numerical  results  in  Fig¬ 
ures  7,  9  i  l1*  and  15. 

Of  course  further  experimental  work  must  1  i  done  to  establish  the 
practical  use  of  this  technique,  especially  as  regards  below-surface  cracks 
that  do  not  penetrate  the  solid.  The  results  do  indicate  the  potential 
for  such  a  technique.  It  is  a  visual  method  whose  features  change  qual¬ 
itatively  as  well  as  quantitatively  when  a  crack  Interupts  the  flow  of 
induced  currents. 
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Stream  lines  for  induced  currents  in  a  cracked  square  plate 
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Figure  8.  Stream  lines  for  induced  currents  in  a  cracked  square  plate 


bi  H 


Figure  10.  Self  induced  field  et  the  center  of  the  induction  coil  u  e 
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function  of  coil  poeition.  a0  *  4,  Xj  ■  0,  ho  »  1. 
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Figure  12.  Eddy  current  intensity  factor  at  crack  tips  as  functions  of 
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coil  position,  aQ  «  4,  ■  0,  hQ  *  1. 
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